MAX-PLUS DEFINITE MATRIX CLOSURES AND THEIR 

EIGENSPACES 



SERGEI SERGEEV 

Abstract. In this paper we introduce the definite closure operation for matri- 
ces with finite permanent, reveal inner structures of definite eigenspaces, and 
establish some facts about Hilbert distances between these inner structures 
and the boundary of the definite eigenspace. 



1. Introduction 

This paper is a contribution to geometrical understanding of some algebraic re- 
sults on max-plus eigenspaces that were obtained by P. Butkovic in (see also [3] ) . 
The sources of geometrical inspiration for our work are [5] and J3j , as well as |S] 
and UJ. Our approach is closer to that of [8] and since we think of max- 
plus semiring with its simplifying total order rather than of generalized algebraic 
structures of idempotent analysis. However, our viewpoint differs from that of jHj 
and [2j in that we use basic tools of max-algebra instead of the more sophisticated 
machinery of convex geometry. 

The paper is organized as follows. In Sect. we recall the basic tools of max- 
algebra that we need. In Sect. we introduce definite forms for max-plus matrices 
with nonzero (finite) permanent and prove that closures of all definite forms of a 
given matrix coincide. Thus we can introduce the 'definite closure' of any max-plus 
matrix with nonzero permanent. We also introduce definite eigenspaces, make some 
observations on systems of inequalities that define them, and consider an application 
to the cellular decomposition introduced in In Sect. ^ we use a representation, 
due to 2 , of the definite eigenspace that reveals some inner structures. Then we 
establish some facts about Hilbert distances between these inner structures and the 
boundary of the eigenspace. 

The author wishes to thank P. Butkovic, A. Churkin, A. Kurnosov, A. Sobolevskh, 
and anonymous referees for helpful comments on the paper. 

2. Some tools of max-algebra 

In its ordinary setting, max-algebra is linear algebra over the semiring R max - This 
semiring is the set R U {— oo} equipped with the operations of 'addition' © = max 
and 'multiplication' = +. Its 'zero' is equal to — oo and its 'unity' 1 is equal 
to 0. 

The semiring R max resembles R + , the positive part of the field of real numbers: 
in both structures multiplication admits inverses and enjoys distributivity over 
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addition, but subtraction is not allowed. In R max , however, the addition is 
idempotent, i.e. for any a € R m a X we have a ffi a = a. In R + this is certainly not 
the case. The semiring R max can be seen as an idempotent 'dequantization' of K. + 
(see, e.g., QJ] and [TT]L 

In max-algebra, it is also possible to exponentiate and to take roots. These 
operations are nothing but conventional multiplication and conventional division, 
respectively. Indeed, for any a / one has a n = a x n and a~ = a/n (for the 
remaining case n = 0™ = since we assume that a = for any a) 

One of the principal objects of max-algebra is R„ lax , the set of n-componcnt 
vectors with components in R max . This set is equipped with 'addition' (a; © y)i = 
Xi © yi and with 'multiplication' by any max-plus scalar (i.e. by any element of 
E max ) (a y)i = a tji. The set R" ax equipped with these operations, as well 
as subsets of M^ ax closed under these operations, will be called max-plus spaces. 
Below the notation will be frequently omitted. 

These max-plus spaces resemble linear spaces in that the laws of associativity 
and distributivity hold, but again there is no subtraction and there is idempotency 
of addition. Structures of this kind are called idempotent semimodules and are a 
central object of the study in the idempotent analysis, see 

A max-plus space S C KSiax i s sa, id to be finitely generated if there is a set of 
vectors { Ui, . . . ,v s } such that for any y 6 S one can find scalars a\, . . . , a s such 
that y = @J =1 onvi. The set { v\, . . . , v s } is the generating set of S. It is minimal 
if no Vi can be expressed as a linear combination of the other generators, i.e. if 
there are no equalities of the form 

(!) ».=0¥j' 

The minimal generating sets will be also called bases. 

The following crucial result is due to Moller and to Wagneur |17l IT%| (it is 
also contained in [H] and |10j)- 

Proposition 1. // {u\, . . . , u s } and {vi, . . . , Vt} are two bases of a max-plus space, 
then s = t and there is a permutation a and a set of nonzero scalars {oti, . . . , a s } 
such that Ui — aiV a ^ for all i = 1, . . . , s. 

Prop, n means that if we have found a finite base for a max-plus space, then it 
is in some sense unique: we can only multiply the vectors of the base by nonzero 
scalars. For more information on max-plus bases we refer the reader to |7j. 

The max-plus matrix algebra is formally analogous to the conventional matrix 
algebra (minus subtraction and plus idempotency): (^4 © B)ij — Aij © Bij and 
(AQB)ij =© fc 4 0% 

Let us introduce two important characteristics of max-plus matrices. The first 
characteristic deals with the cyclic permutations. Let A be an n x n max-plus 
matrix. Here and below N will stand for the n-element set {1, . . . , n}. Denote by 
C n the set of all cyclic permutations r that act on the subsets of the set N. For 
t G C n denote by K(t) the subset on which r acts and by | K(t) \ the number of 
elements in this subset. Then 
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is the maximal cycle mean of the matrix A (the notation A (A) for the maximal cycle 
mean of A will be used throughout the paper). The summand (Qi e xM ^ir(i)) |ir(T) 
is called the cj/de mean of the cyclic permutation r. 

The cyclic permutation whose cycle mean is maximal will be called critical. 

The second characteristic deals with the permutations of N. For any square 
n x n max-plus matrix A, its permanent is defined as 

n 

(3) per(A) = O^-W- 

crgS„ i=l 

Here 5„ is the group of all permutations of N. 

The summand 0"=i -^-ia-U) is called the weight of the permutation a, so the max- 
plus permanent is the maximal weight of all permutations of N. The permanent of 
A is said to be strong (following 3 ), if A has only one maximal permutation, i.e., 
only one permutation with maximal weight. 

In the papers [8] and [5] the max-plus permanent is called the 'tropical deter- 
minant'. To some extent the max-plus permanent can overtake the role that the 
usual determinant plays, see |3] and the just mentioned papers for details. We also 
refer the reader to pQ and |10j . for a symmetrized version of max-algebra, which 
admits subtraction and determinants. 

An n x n max-plus matrix A is called invertible, if there is another n x n max- 
plus matrix B such that the products AB and BA are both equal to the max-plus 
identity matrix /. Prop. ^ implies that A is invertible iff there is a permutation a 
and a set of nonzero scalars a±, . . . , a n such that 

A = \a> t , ifj = o-(*); 

lo. otherwise. 

So the class of invertible matrices in max-algebra is very small. But it makes 
sense to calculate the series 

(4) A* = I © A © A 2 © . . . , 

where I is the max-plus identity matrix. If the sum of this series exists, then it is 
called the max-plus algebraic closure of A. It is an obvious analogue of (I — A)^ 1 . 

Powers of max-plus matrices and max-plus closures play important role in opti- 
mization on graphs. Indeed, if we associate an n-node graph with an n x n matrix 
A and let Aij be the weight of the edge of this graph, then the entry (A m )ij 
of the matrix A m will represent the maximal weight of paths of length m running 
from i to j. The entry (A m )u is the maximal weight of all cyclic paths, i.e. cycles, 
of length m that traverse i. In other words, it is the maximal weight of the cyclic 
permutations r such that i € K(r) and | K{t) |= m. Analogously, (A*)ij, for 
i =/= j, is the maximal weight of all paths of any length running from i to j. The 
path from i to j whose weight equals {A*)ij is called optimal. 

The following proposition, due to Carre j^, solves the problem of existence of 
the closure. 

Proposition 2. The closure of the matrix A exists if and only if X(A) ^ 1. 
Max-plus closures enjoy the property 

(5) (A*) 2 = A*. 

Otherwise stated, the inequality A*jA* k < A* k holds for all i,j, and k. 
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As A* is an equivalent of (I — A) 1 , some algorithms of linear algebra can be 
adapted to calculate closures (that is, all optimal paths on a graph) in max-algebra, 
see 01, H3|, and H5|. 

Another important tool is the max-plus spectral theory. For A, an n x n max- 
plus square matrix, the max-plus spectral problem consists in finding an n-element 
vector x, such that not all of its components are 0, and a scalar A such that 

(6) Ax = Xx. 

The scalar A is a max-plus eigenvalue, and the vector £ is a max-plus eigenvector. 
The set of max-plus eigenvectors associated with A is closed under addition and 
multiplication by any nonzero scalar. Therefore it is called the max-plus eigenspace 
(associated with A). The eigenspace associated with the maximal eigenvalue of A 
will be denoted by eig(A), and the space generated by the columns of A will be 
denoted by span(A). 

One of the main results on max-plus spectral theory is the following. 

Proposition 3. The maximal eigenvalue of any max-plus square matrix is equal 
to its maximal cycle mean. 

The eigenspace associated with the maximal eigenvalue is easy to describe. 

Proposition 4. Let A be an n x n square max-plus matrix such that X(A) = 1. 
Then 

1) eig(A) is generated by the columns A*- of A* such that AA*^ — A*^; 

2) AA*j = A*j iff there is a critical cyclic permutation t such that j s K(t). 

Columns of A* corresponding to vertices of the same cycle are proportional (in 
the max-plus sense) to each other. 

Proposition 5. Let A be an n x n max-plus matrix such that X(A) = 1 and let r 
be a critical cyclic permutation. Then for any I 6 N and i € K(t) we have 



If the graph associated with A is strongly connected (i.e. for any i and j there 
is a path with nonzero weight running from i to j), then A is said to be irreducible. 
In this case the maximal cycle mean of A is known to be its only eigenvalue. In the 
reducible case this eigenvalue need not be unique. If another eigenvalue exists, then 
it is the (only) eigenvalue of some maximal irreducible submatrix of A. However, 
not all maximal irreducible submatrices of A yield eigenvalues for A. 

For more details on max-plus spectral theory, as well as for the proofs of the 
above propositions, we refer the reader to pQ, [H], ^H], and [T§| . 



We begin this section with the following important proposition. The proof makes 
use of the uniqueness of the base (Prop.^l. 

Proposition 6. Let A and B be square max-plus matrices such that X(A) ^ 1, 
and X(B) ^ 1. Then span{A*) — span(B*) if and only if A* = B* . 



(7) 
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Proof. First let us prove that whenever 

(8) a* 0".i> 

iiti 

the column A*± is proportional to A*j for some j =/= i. 

Suppose JSJl holds. Then there is an I such that A* { — 1 = and that 

A* u > aj. Combining this we obtain that A* a A* H ^ 1, hence A* a A* u = 1 (otherwise 
there is a cycle whose weight exceeds 1). This means that there is a critical cycle 
with weight 1 traversing i and I, and due to Prop. A* t — aiA* r We conclude 
that no column of A* can be expressed as linear combination (JHJ without being 
proportional to some of the columns involved in this combination. 

Further let {A* ri , . . . ,A* rk } be the base of span(A*) = span(£T). If we use 
columns of B* to form the base, then, due to Prop. it must be of the form 
{ B* , . . . , B* j, so that B* = aiA* r for i — 1, .... ft and some nonzero aj. 
All remaining columns of A* (of B*) are proportional to base columns of A* (of 
B*). So every column of A* is proportional to some column of B*, and vice versa. 
This implies that the rows A*, and A*, are proportional iff the rows B* and B*. are 
proportional. 

Let the columns A\ and A*- be proportional. Then A*jA* { = 1, hence there 
is a critical cycle containing i and j. Due to Prop. [S] the rows A*, and A*, are 
proportional, so are the rows B* and B*. , and, again due to Prop. 03 so are the 
columns -B* and £?* . We conclude that the columns A* and A*j are proportional 
iff so are the columns B* and B%. 

Now it is clear that { A* ri , . . . , A* Tk } is the base of span(A*) iff { B* ri , . . . , B% k } 
is also the base, so that B* r . = onA* ra . We can assume w.l.o.g. that r; = i. 
Consider the decomposition of a into cyclic permutations and let r be one of them. 
Then B*,^ = a % and A* T ^ = a^ 1 for all i e K(r). If O ieX (r) a i > tnen 5 * nas 
a cycle with weight greater than 1, and if O ieA: ( r ) a-i < 1, then so does A* . The 
only remaining possibility Oieif(r) ai = ^ implies that all columns of A* and -B* 
with indices belonging to K(t) are proportional. Then K(t) must be a singleton 
for any r, otherwise the minimality of the bases is violated. This implies that a is 
the identity permutation, and A* = _B* for any i = 1, . . . , k. Taking into account 
that there is a one-to-one correspondence between the sets of proportional columns 
of A* and B* , and that all columns of A* and B* are proportional to some base 
columns, we conclude that A* = B* . □ 

Now consider matrices with maximal cycle mean equal to 1 and with all diagonal 
entries equal to 1. Following [3], we call such matrices definite. The following 
proposition contains some simple facts on eigenspaces of such matrices. The third 
statement is an easy consequence of , Ch. IV, Th. 2.2.4, and its proof is recalled 
here for convenience of the reader. We consider the general reducible case, in which 
the eigenvectors may have zero entries. For any y S K" ax , the index set K such 
that yi 7^ iff i G K is called the support of y and is denoted by supp(y). 

Proposition 7. If A is a definite matrix, then 

1) it has a unique eigenvalue equal to 1; 

2) eig(A) = span(A*); 

3) an eigenvector with the support K C N exists iff Aij = for all i € N\K 
and j G K . 
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Proof. 1) Any eigenvalue of A is the maximal cycle mean of some of its subma- 
trices, and the maximal cycle mean of any submatrix of A is equal to 1. 

2) follows from Prop.^J 

3) If an eigenvector x such that supp(rr) = K exists, then A^jXj = for all 
i G N\K and all j £ K. Hence Aij = for all such i and j. 

Conversely, if the set K satisfies the condition, we look for eigenvectors x such 
that supp(.x) = K . We may reduce the system Ax = x to the system Akk%k — %K> 
where Akk is a submatrix of A standing on the rows and columns with indices 
belonging to K, and xk is a vector with | K | nonzero components. The space 
eig(^4 kk) is generated by the columns (Akk)^, where j £ K. Taking any com- 
bination of all these generators with all coefficients not equal to we obtain an 
eigenvector with the support K. □ 

In the case of definite matrices we have one more implication of the uniqueness 
of the base. The proof is similar to that of Prop. 

Proposition 8. If A is definite and span(A) = eig(A), then A = A* . 

Proof. Let { A. S1 , . . -A. Sk } be the base of span(A) = eig(A) = sp&n(A*). Due 
to Prop, n if we use columns of A* to form the base, then it must be of the form 
{ A* tl , . . . A* tk } so that A, Si = otiA* t . for i = 1, . . . , k and some nonzero on. More 
precisely, a, = A tiSi , and this implies A tiSi A*. t . = 1. Then there is a critical cycle 
traversing s, and U, so A* s . — A tiSi A* t . according to Prop. El So A. Si = A* s . for all 
i = 1, . . . , k. 

Now assume that there are columns A.j and scalars a.i such that 



and that no column of the base is proportional to A.j. It follows from @ that 
there is an index m such that Ajj = 1 = a rn Aj Sm . The columns A.j and A. Sm 
are not proportional, hence there is an I such that Aij > a m Ai Sm . This implies 
AijAj Sm > Ai Sm and A^ s > Ai Sm . This is a contradiction, since we have proved 
that A* s . = A. Si for any i = 1, . . . , k. So any column of A is proportional to some 
column of the base. But if A.i and A.j are proportional, then AijAji = 1, hence 
A*i and A*j are also proportional with the same coefficient. This implies A = A*. 



Let us now introduce a definite form of a matrix. Consider an n x n max- 
plus matrix A that has nonzero permanent, i.e., at least one permutation whose 
weight is not equal to 0. Let a be one of the maximal permutations of A. If this 
permutation is not the identity permutation, then we turn this permutation into 
the identity permutation by rearranging the columns of A. Then we divide (in the 
max-plus sense) all columns by the corresponding diagonal entries, thus obtaining 
a matrix A' with entries 



k 



(9) 



A.j — ct{A. Si , 



□ 



(10) 



-1 



Obviously, passing to a definite form of a matrix does not alter its span: span(A') = 
span(A). 
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It is clear that A' is definite, hence we call it the definite form of A corresponding 
to the permutation a. For example, if 




then 

/ -oo -5^ 
A' = 1-4 -4 
\-oo -5 

is the definite form of A corresponding to u = (231). (Here and in the sequel we 
prefer not to use the symbols and 1 in numerical examples.) 

However, in general there are many maximal permutations and many definite 
forms corresponding to them. The fact that we want to prove is that closures of all 
definite forms coincide. It is convenient to pose the problem as follows. Let A be 
definite, let it have maximal permutations different from the identity permutation, 
and let a be one of them. Let A' be the definite form of A corresponding to a, 
then its entries are defined according to (|10(l . A' also has maximal permutations 
different from the identity permutation, and cr _1 is one of them, since 

(11) ^-Hi) = K-m)i- 

First we prove the following proposition. 

Proposition 9. Let A be a definite matrix and let a be one of its maximal permu- 
tations. For A' the definite form of A corresponding to a, eig(A') = eig(A). 

Proof. Consider the decomposition of a into cyclic permutations. Let r be 
the number of these permutations and denote these permutations by t\ (where 
I = 1, . . . ,r). Denote by K(t{) the index set on which 77 acts. The sets K{t{) are 
pairwise disjoint and |Jj=i K( T i) = N ■ 

We prove the inclusion eig(A) C eig(A'). Let y be an eigenvector of A with the 
support M C N. Due to the third statement of Prop. the set M must be of the 
form lj zgi K{ti) for some L C {1, . . . , r}, and we have Aij = for all i 6 N\M 
and j S M. 

Now note that every cyclic permutation 77 is critical, hence for any 77 such that 
Ki G M we have 

(12) O A in(i)yn{i)yi l = 1 

ieK(n) 

On the other hand we have A iTl ^y Tl ^ ^ y^. This inequality is an equality, 
otherwise the violation of (|12f) occurs. Thus we obtain 

(13) Aa(i)ya(i) = Vi 

for all i e M. 

Now note that, since e M for all j S M, we have that A[- — for all 
i G N\M and j G M, the same as for Aij. Therefore it suffices to prove that 
if AmmVm = Vm, then A' mm dm = Dm- In other words, we want to show that 
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A'^yj < y l for all i,j E M. But due to (JTHJ) and lfT3|) 

A 'ijVj = A i<y{j) A ja(j)Vj = 

and the proof of the inclusion eig(A) C eig(A') is complete. The opposite inclusion 
is proved analogously, after passing to the definite form associated with a (the 
matrices A and A' will interchange). □ 
The fact that we want to prove is now an immediate consequence of Proposi- 
tions HO and El 

Proposition 10. Closures of all definite forms of any matrix with nonzero perma- 
nent coincide. 

We have just proved that closures of all definite forms coincide. Now we can 
define the definite closure of any nxn max-plus matrix with nonzero permanent to 
be the closure of any of its definite forms. Due to the second statement of Prop.0 
for any definite form A' of A we have eig(A') = eig((A')*) = span((A')*). So the 
eigenspace of the definite closure of A coincides with the eigenspace of any of its 
definite forms (they are all the same), and is generated by columns of the definite 
closure. The eigenspace of definite closure will be called the definite eigenspace. 

The third statement of Prop.[7]suggests that, if we want to work with the space of 
eigenvectors of definite closure of non-full support, then we must confine ourselves 
to the corresponding submatrix. 

Further we always assume that the eigenvectors considered have full support, 

i. e., that we study the eigenvectors with certain support and have passed to the 
corresponding submatrix. 

Let us show that the definite eigenspace can be described by some system of 
inequalities. 

Proposition 11. Let Abe a definite matrix. Then its eigenspace (and the eigenspace 
of its closure) is the set X — { x \ A^ ^ x^xj , i ^ j, A^ — A*- } 

Proof. Let x be an eigenvector of A corresponding to its maximal eigenvalue 
1. It satisfies the system Ax = x, therefore it satisfies all inequalities of the form 
Aij ^ XixJ 1 . Hence x 6 X. 

Conversely, if all inequalities Aij ^ XixJ 1 , for i ^ j and Aij — A*j, are satisfied, 
then absolutely all inequalities A^ ^ satisfied. Indeed, if A^ < A*^, then 

the optimal path from i to j is not the edge but it traverses other nodes, say, 

ii, . . . , i k . Then A* 3 = A Ul . . . A lkj where A ih = A* n , . . . , A ikJ = A* kJ (all edges 
(i,ii), . • • , (ik,j) are optimal paths). The inequality A^ ^ XixJ 1 is now an easy 
consequence of the inequalities Au 1 ^ XiX^ , . . . , Ai k j ^ Xi k xJ 1 that are satisfied. 
So all inequalities A^ ^ satisfied and this implies Ax — x. □ 

An inverse problem can also be posed. Suppose that there is a system of inequal- 
ities {aij ^ XixJ 1 }, with at most one inequality per each pair (i, j), and the set of 
vectors with full support defined by this system is not empty. Is it a full-support 
subspace of an eigenspace of a definite matrix? 
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To answer this question, consider the matrix A whose entries Aij are equal cither 
to 1, if i = j, or to a^, if there is an inequality of the form Xixj 1 , or to 0, if 

there is no such inequality. Then we have the following proposition. 

Proposition 12. The set X = {x \ A%j SC XiXj 1 } is nonempty if and only if A is 
definite. 

Proof. If A is definite then the set of its eigenvectors associated with the 
eigenvalue 1 is nonempty and, due to Prop. 1111 it is precisely the set X (some of 
inequalities being redundant). 

Conversely, let X be nonempty, then there exists x G X. Take an arbitrary 
cyclic permutation r. As a consequence of all inequalities A iT ^ x i X T(i)' wnere 
i £ K(t), we obtain that OiejfYr) -^«t(*) ^ Hence all cycle means of A are not 
greater than 1, and A is definite. □ 

Now consider the following application. 

Let V be an m x n max-plus matrix with at least one nonzero entry in each 
column, and let y £ R™ H have full support. Then, following we can define the 
combinatorial type of y with respect to V. It is an m-tuple S of subsets Si, ... , S m 
of {1, ... , n}, such that i S Sj whenever VkiVf^ — VjiVj 1 - It is proved in |U] 

that the collection of the sets 

(14) X s = { y | V ki V^ < ykyj 1 for j, k = 1, . . . , m and i e S,- } 

defines a cellular decomposition of the full-support subspace of K™ ax - 
Consider an m x m matrix l^" 5 whose columns are defined by 



ieSj 



(15) V.f = ^K,, 
if Sj is not empty, and by 

(16) Vg 




if is empty. Then it is clear that 

(17) X s = {y\V? ^ yi yr x for i, j = l,...,m}, 

and we immediately have the following proposition. 

Proposition 13. 1) The cell Xs is nonempty if and only if X(V S ) ^ 1; 

2) IfX(V s ) < 1, then I s = {(/£ span((V s )*) | supp{y) — N }; 

3) // \{V S ) < 1, tAen X s = { x \ V? < Xi xj\ % ± j, V? = (V% }. 

We close this section with three examples. 



Example 1 Consider the matrix 
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The only maximal permutation of A is (231), and the only definite form is 

/ 1 N 

A' = -3 2 

\-4 -5 ; 

The definite closure of A is 

/ 2^ 

(A')* = -2 2 

\-4 -4 0y 

Fig. 1 displays the cross section by z — of span(A) (left) and span((A')* 
eig(A') (right). 




4 x 



Figure 1 . The definite closure operation for A 

Note that eig(A') is the set 

{(x, y, z) | xy' 1 ^ 0, yz" 1 ^ 2, zx" 1 ^ -4 }, 
in accordance with Prop. ITTI 



Example 2 Consider the matrix 

I' 1 (1 




It has two maximal permutations: (13) (2) and (231). Therefore it has two 
definite forms, namely 



and 




But, in accordance with Prop. El the definite closure 

/ -1 2 N 

(B'y = {B"y = i 03 

1-3 -4 0/ 
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is unique. 

Fig. 2 displays the cross section by z — of span(S) (left) and span((_B')*) 
eig{B') (right). 



4 x 



4 x 



Figure 2. The definite closure operation for B 

In accordance with Prop. ^2 the space eig(S') is the set 

{ (x, y, z) | yx- 1 > 1, xy' 1 > -1, xz~ x > 2, zy' 1 > -4 }, 

or, equivalently, the set 

{ (x, y, z) | yx~ x ^ 1, xy^ 1 > — 1, xz^ 1 ^ 2, zx~ x ^ -3 }. 

The first system of inequalities corresponds to the definite form B' and the second 
one corresponds to B". 

Comparing Fig. 1 with Fig. 2 we see that eig(A') has 'interior' whereas eig(-B') 
does not have 'interior' (for the exact meaning of the term 'interior' see Sect. 01 
below). As a consequence of 0, Th. 4.2, or [H], Th. 4.2, one can obtain that 
eig(A'), for A' a definite form of A, has 'interior' if and only if A (or equivalently 
A') has strong permanent. This fact will be revisited in Prop. H*4l of this paper. 



Example 3 Consider the matrix 

/l 4 6 7\ 
7=4 1 5 8 . 
\0 0/ 

Let S = ({2,3},{4},{1}), P = ( 0, {4}, {1, 2, 3} ), U = ( {3}, {1, 3, 4}, {2} ), and 
W = ( {1, 4}, {2}, {3} ) be four combinatorial types. Do they exist in the cellular 
decomposition? If they do, what vectors generate the respective cells Xs, Xp, Xjj, 
and Xyjl 
For S: 

/ -1 1\ / -1 3 

V s = -1 4 1, (V s )* = 4 
\-4 -8 0/ \-4 -5 

Hence X s exists and is generated by [4 4 0] T , [4 5 0] T , and [3 4 0] T . 
For P: 

/ -1 6\ 

V p = -co 5 = (V p )*. 
V-oo -8 0/ 
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Hence Xp is generated by [0 — oo — oo ] T , [7 8 0] T and [6 5 0] T . However, 
[0 — oo — oo] T does not have full support and does not belong to Xp. 
For U: 

( 1 4\ / 1 4\ 

v u = -l o 1 , (v u y = -l o 3 . 

\-6 -4 0/ \-5 -4 0/ 

Hence X' 7 exists and is generated by [5 4 0] T and [4 3 0] T . 
For W: 

/ 3 6 
V w — | 3 5 
\-l -1 

We have X(V W ) — 3 > 1, hence W does not exist in the cellular decomposition. 

Fig. 3 displays span(V) (blue), X$ (red), Xp (light grey), and X;y(dark green), 
projected onto z = 0. The generators of span(V) (the larger circles) and the 
generators of X$ and Xjj{the smaller squares) are also shown. 




2 4 6 8 x 



FIGURE 3. Three cells of a cellular decomposition 



4. Inner structures of definite eigenspaces 

We have introduced the definite closure operation and have given an external 
description of a definite cigenspace in terms of a system of inequalities. In this 
section we give an internal description of definite eigenspace and measure Hilbert 
distances between the structures involved in this description and the boundary. 

For further considerations we need the following notions and notation. Let A be 
a definite matrix. The sets Xij — {x \ Aij = Xixj 1 } will be called the supporting 
planes of eig(A), and the sets = PI eig(A) will be called the faces of cig(A). 
The boundary, i.e., the union of all faces of eig(^l) will be denoted by T(eig(A)). 
The set of eigenvectors not belonging to the boundary will be called the interior of 
eig(A) and denoted by int(eig(A)). 

Also, denote by A the matrix obtained from A by replacing the diagonal l's by 
O's, and denote by A^ the matrix obtained from A by replacing the diagonal 
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O's by l's. For example, if 



then 



and, e.g., 

K -4 -6 

Also, the maximal cycle mean of A is \(A) = —1.5, and 

/ -2.5 2.d> 
A x = 2.5 
\-3.5 -5.5 

(we write A\ instead of ^a(A) f° r * ne sa ^ e °f simplicity) 
The following crucial proposition can be derived from 

Proposition 14. Let A be a definite matrix such that X(A) ^ 0. 

1) If A does not have strong permanent then eig(A) does not have interior; 

2) If A has strong permanent then eig(A) has interior, and 

(18) mt{eig(A))= (J eig(A^). 

A(A)^A!<J 

Proof. 1) If A does not have strong permanent, then it has a maximal permu- 
tation that differs from the identity permutation. Let a be such permutation, and 
assume that there is a y belonging to int(eig(^4)). Then A ia ^y a ^ < yi for all i. 
After multiplying all these inequalities and then cancelling the product Q - yi one 
obtains Oi^-io-(i) < 1- This implies that a is not maximal, a contradiction. 

2) y belongs to int(eig(A)) if and only if J y i ^4ij%- < y% for all i. This takes 
place if and only if there is a [i < 1 such that AijUj ^ W«i or i equivalently, 

A^y — y. If we take /i A(A), then 1 is the maximal cycle mean of A^, hence 
eig(^4 M ) exists. The proof is complete. 

□ 

If X(A) — 0, i.e. if the graph associated with A is acyclic, then representation 118(1 
is replaced by the representation 

(19) mt(eig(A))= (J eig(^). 

0<^i<1 

In the sequel, we always assume that A is definite and has at least one off-diagonal 
entry not equal to 0. 

According to Prop. 1111 the eigenspace eig(A) is the set 

x = { x I A ij < XixJ 1 , i ^ j, A i0 = A*j }. 

Analogously, the eigenspace eig(A fl ), for any ^ involved in (|18|1 or I)1H[I. is the 

set 
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X H = { x I M ^ XiXj \ i^j,fl 1 A lJ = (A*)ij }. 

We need the following proposition mainly for the proof of Prop. El 

Proposition 15. Let n be a scalar such that \(A) ^ fi < 1, if A(A) > 0, or such 
that <»<!, if A(A) = 0. If (A*) y = (A M ) tf , ften Ay = Ay. 

Proof. In both cases considered the maximal cycle mean of A M is equal to 
1, hence A* exists. Let A*j > Aij, then A*» = Au x ...Ai k j for some ii,...,ik 
not equal to i or j. Since /i < 1, we have fj, Aa t . . .fx Ai k j > n~~ 1 Aij, hence 
(A*)^ > (A M ) y , ' ' " □ 

The eigenspaces eig(A jU ) are the inner structures mentioned above. Now we 
are going to measure the Hilbcrt distances between these inner structures and the 
boundary T(eig(A)). 

The Hilbert distance between the two vectors x and y both having support K is 
defined to be 

(20) d H (x,y)= x l xj 1 y^ 1 y J . 

i,jeK 

Note that in [S] the Hilbert distance is defined as an inverse of the quantity 
d}j{x,y)- If the supports of x and y differ, then we assume the Hilbert distance 
between x and y to be infinite. 

It can be easily verified (see also jS], Th. 17) that the following properties hold: 

1) dn{x, y) }z 1, and du(x, y) = 1 iff x = \y, where A is a nonzero scalar; 

2) d H (x,y) = d H (y,x); 

3) d H (x,y)d H (y,z) ^ d H (x,z). 

In fact these properties show that dn is a semidistance (recall that 1 = and 
= +). Indeed, dn{x, y) = 1 = whenever x is equivalent to y modulo 

x ~ y 3A ^ : x = Xy. 

This semidistance is induced by the range seminorm 

(21) \\x\\ = x iX j\ 

introduced in see also [7j- However, by a slight abuse of language we will refer 
to dn as to a distance. 

Now we measure the distance between an arbitrary y £E eig(A) and the supporting 
plane X^ = { = Aij }, i.e., the minimal distance between y and x S X^. 

From 5 , Th. 18 it follows that this minimum is attained at the maximal vector of 
Xij not greater than y. Denote this vector by y l K Its coordinates are very easy to 
find: 

Vi =(Bi x i\xi^yi} = yi for I ^ i, j; 

(22) yl 1 =(B{xi\xi^ yuA^Xi < y } } = Ayy,-; 
Vj = ©{ Xj I Xj < yj,AijXj ^ yi } = yj. 

The distance ll'iPp between y and X^ is then equal to 

(23) d H (y,X ij )=AT j 1 yT 1 y i . 
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However, what we need is the distance between y and r(eig(A)), i.e, the minimal 
distance between y and Tij. The following proposition makes our life simpler. 

Proposition 16. The distance between y £ eig(A) and the boundary T(eig(A)) is 
equal to the minimal distance between y and supporting planes. 

Proof. Clearly the minimal distance between y and supporting planes is not 
greater than the distance between y and T(eig(^4)). Suppose Xy is the supporting 
plane such that the distance between y and X^ is minimal. This distance is equal 
to the distance between y and y % K If y l i belongs to eig{A) and hence to then 
we are done. Suppose not; then the system of equalities Akiy] = y k must be 
violated for some k £ N. Note that, if k ^ i, then yV = yk (see 11220 ). and there is 
no violation. So the violation must take place for k — i. There must be an / such 
that Auyf > y\° , i.e. such that Auyi > A^yj. Now consider z such that Zi — Auyi 
and Zk = yu for any k ^ i. Then z belongs to and dn{y, z) — A^ y~ y i4 This 
distance is strictly less than the distance between y and , a contradiction. □ 

Consequently, 

(24) d H (y,T(eig{A))) = f\ A^yfy,. 

The key idea of Prop.lT7lbelow is that X(A)~ 1 , if X(A) is invertible, is the largest 
radius of Hilbert balls contained in eig(A). It can be said that A(^4) _1 is the radius 
of inscribed Hilbert balls, as depicted on Fig. 5. 

Let r be any critical cyclic permutation of A. 

Proposition 17. 1) In the case X(A) > for any y £ eig(A) the distance 
between y and T{eig{A)) is not greater than X(A) ; 

2) Let [i be such that X(A) ^ ji < 1, if X(A) > 0, or such that < [i < 1, if 
X(A) = 0. Then for any y £ eig(A fJj ) the distance between y and T(eig(A)) 
is not less than fi ; 

3) In the case X(A) > 0, for any i,j £ K(t) such that j = r(i), and any 
y £ eig{A\), the distance between y and the face is equal to A(A) _1 . 

Proof. 1) The distance between y and T(eig(A)) does not exceed the minimal 
distance between y and supporting planes Xy that correspond to the edges (i, j) 
of the cyclic path determined by r, and this minimal distance is not greater than 
X(A)-h 

d H (y,T(eig(A))) < /\ ieK(r) Ar^y^ < 

^(O i eA' M A> 1 W )^=A(i)- 1 . 
2) If y £ eig(A M ) then, since eig(A l _ l ) = span(A*), we have 

(25) y = a k (A;). k , 

where M = supp(a). Substituting 125|l into (|24Jl . we get 

(26) d H (y,T(eig(A))) = f\ Ar. 1 /\ ^(A^ 

i^j^Aij^O keMj I 



16 



SERGEI SERGEEV 



Here by Mj we denote the set Mnsupp((A*)j.). Now we estimate (12 fill from below 
and use the inequalities Ay ^ /x(A*)y and (A*)ij(A*)jk ^ (A*)ifc (see ©): 

d H (y,r(ei g (A))) ^ A,.-,..» «,A*, M ^(^Wtt > 

55 A,.-,.» ...A^.u ^(^'WjiWa ^ 

3) The distance between y G cig(AA) = span(A^) and the supporting plane 
is equal to 

(27) d H (y,X l3 ) = f\ al\Al)^ o^aV 

keMj leMi 

The cyclic permutation r of A\ has the weight 1. Hence for all i,j G K(t) 
such that j = r(i) we have, according to Prop. 03 that Aij = A(A)(At)jj and 
(Al) ij (Al) jl = (A* x )u. Note that {A* x ) tj ^ for all i,j G if and therefore (A^) j7 = 
if and only if (A^)^ = 0, i.e. Mj and Mj coincide. Making use of all this we 
write the upper estimate for du(y, X^j): 

= ® ieMi HAyHAi^HADj'iAiu = x(Ay\ 

We also have dn(y, r(eig(A))) ^ X(A)^ 1 and therefore (see Prop. lPo|) dn(y, Xy) = 
d H (y,T ij ) = X(A)- 1 . □ 

The sets r(eig(A M )), for fi < 1, are the subsets of eig(A) equidistant from 
r(eig(A)), as Prop. IT%) suggests. 

Proposition 18. For all /i such that X(A) ^ /i < 1, if X(A) > 0, or such that 
< /i < 1, if A(A) = 0, the distance dn{y, T(eig(A))) is equal to [iT x if and only 
tfyeT(eig(A^)). 

Proof. If X(A) > and [i = A(A) then the statement readily follows from the 
observation that A\ does not have strong permanent and therefore (see Prop. I14|) 
eig(A>) does not have interior. 

Let us consider /j, > X(A). First, the equality 

AijVjVi 1 = A 1 

implies A^y — y. So, if dn (y, r(eig(A))) = n^ 1 then y G eig(A M ). Assume that y 
belongs to the interior of eig(A M ). Since A^ is definite and has strong permanent, 
we can use representation (I18|) or 119fl and obtain n < 1 such that y G eig(A |llft ). 
Now statement 2) of Prop. IT7I implies that dn(y, r(eig(A))) ^ (/ik) _1 > [i. This is 
a contradiction, so y G r(eig(A p )). 

Suppose now that y G r(eig(A M )). It means that there are i ^ j such that 
ViVj 1 — (Afj,)ij, where (A*)^ = (A^y. According to Prop. El this face corre- 
sponds to the face of eig(A) determined by the entry Ay, and the distance between 
these two faces is clearly n~ x . □ 

Throughout this section we dealt with eigenvectors having full support. But let 
us recall the third statement of Prop. [7| It says that there might be eigenvectors 
with nontrivial support K. The distance between these eigenvectors and part of any 
face with full support would be infinite. Also, these eigenvectors are eigenvectors of 
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the submatrix Akk- Therefore it is presumable, in this case, to pose the problem 
of finding d H {y,T(eig(A K K)))- 

We conclude this section with two examples. 



Example 1 In the beginning of this section we considered the definite matrix 







-4 




I 









-7 



with the maximal cycle mean of A equal to A = —1.5. Now we pick the following 
three members of the {A^} family: 



and 



A. 



0.5 






-3.5 1.5N 






1.5 


1.5 




-l — 


-4.5 


-6.5 ) 








( ° 


-2.5 


2.5 


A x 


= 2.5 





2.5 




V— 3.5 


-5.5 













5 x 



5 x 



Figure 4. The space span(A) and the eigenspaces eig(A^) 



The left-hand side of Fig. 4 displays span(A), and the right-hand side displays 
the sets T(eig(A fJi )), for /i = (dark blue) \i = 0.5 (green), fj, = 1 (brown), and \i = 
A = 1.5 (red). The lines corresponding to larger values of y, are given smaller weight 
(to help distinguishing between different values of /i in black and white printing). 
We see that there is an injection of systems of inequalities describing eig(A M ) into 
the system of inequalities describing eig(A) in accordance with Prop. [TBI 

Fig. 5 displays, together with span(A), two Hilbert balls with the maximal radius 
A -1 = 1.5 inscribed in r(eig(/l)). The balls touch the 'cycle faces' {(x,y,z) G 
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5 x 



5 x 



FIGURE 5. Two Hilbert balls inscribed in T(eig(A)) 

eig(A) | yx -1 = 1 } and { (x, y, z) 6 eig(A) | xy -1 = —4 }, in accordance with the 
third statement of Prop. El 



Example 2 Consider a Hilbert ball with radius d centered at {Ax} (A is any 
nonzero scalar). It is the set 



1,3 



y = {y i0^rV 



d 1 Xj.x A . It is easily verified that 



or, equivalently, 

Y = { y yiyj 1 > d 

Denote by D the matrix with entries dij = d^'xix^ 
D = D* . Then it follows from Prop. ITT1 that the Hilbert ball is the eigenspace of D 
and the columns of this matrix are its generators. The maximal cycle mean of D 
is clearly The eigenspaces eig(Du) where d -1 < fi ^ 1 are Hilbert balls with 
radii (fid) -1 centered at {Ax}, and eig(D d -i) is precisely {Ax}. 

For a three-dimensional example, set x = [5 4 0] T and d = 3. Then 

/ -2 2 
D = -4 1 
1-8 -7 




Figure 6. Hilbert balls as eigenspaces 



Fig. 6 displays the sets r(eig(Z? /J ) for ft = (dark blue), /i = —1 (green), and 
/i = —2 (brown) with the same convention about the weight of lines as in Fig. 4. 
These sets are concentric spheres centered at eig(Z?_3) = {Ax} (the large red circle 
in the center of Fig. 6). 
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